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Abstract 

By incorporating the curvature 0,^ as a free parameter, it has been found that the tension 
between the high redshift cosmic microwave background (CMB) shift parameter R(z*) data and 
the low redshift Type la supernova (SNIa) and baryonic acoustic oscillation (BAO) data from 
the combination of Sloan Digital Sky Surver (SDSS) and Two-Degree Field Galaxy Redshift 
Survey (2dFGRS) can be ameliorated, and both SNIa+BAO and SNIa+BAO+CMB favor that 
the decelerating parameter q{z) shows a rapid variation in the sign at the small redshift. In this 
paper, with the Monte Carlo Markov chain method, we re-examine the evolutionary behavior of 
q(z) using the latest observational data including the Union2 SNIa, BAO, and CMB data (R(z*), 
Ia{z*), z*) from Wilkinson Microwave Anisotropy Probe 7-yr (WMAP7). For the BAO data, 
four different data sets obtained from the Two-Degree Field Galaxy Survey, the combination of 
SDSS and 2dFGRS, the WiggleZ Dark Energy Survey and the Baryonic Oscillation Spectroscopic 
Survey, are used. Except for the spatially flat case constrained by SNIa+ the WiggleZ BAO 
data, both SNIa and other BAO+SNIa favor that the present cosmic acceleration is slowing down 
irrespective of whether the spatial curvature is included or not. However, once the WMAP7 CMB 
data is included, observations favor strongly the Lambda cold dark matter model, a spatially 
flat universe, and a speeding-up of the cosmic acceleration. Therefore, the inclusion of spatial 
curvature seems to has no effect on alleviating the tension between SNIa+BAO and CMB in 
clear contrast to the previous work in the literature. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 



Two independent Type la supernova (SNIa) observation groups discovered firstly in 
1998 that the Universe is undergoing an accelerating expansion [1, 2]. This discovery was 
further confirmed by observations on the cosmic large scale structure [3, 4] and the cosmic 
microwave background radiation [5]. At the same time, the observations also indicated 
that the cosmic phase transition from deceleration to acceleration occurs only in the near 
past. Interestingly, using the Constitution [6] and Union2 [7] SNIa along with the baryonic 
acoustic oscillation (BAO) distance ratio between redshifts z = 0.2 and 0.35 [8] and the 
Chevalier-Polarski-Linder(CPL) parametrization [9] for the equation of state (EOS) of 
dark energy, the authors in Refs. [10, 11] recently found that the acceleration of the cosmic 
expansion might be slowing down. However, once the CMB data is included, their results 
turn out to be supporting that the universe is undergoing an accelerating expansion with 
an increasing acceleration. Thus, there appears some tension between the low redshift data 
(SNIa+BAO) and the high redshift (CMB) one. This tension can be reduced by various 
methods, which yield different results on the cosmic expanding history [11]. Through 
combining different SNIa and BAO data, Gong et al [12] found that the systematics in the 
data sets does influence the fitting results and leads to different evolutional behavior of the 
decelerating parameter. The effect of different parametrizations of the EOS of dark energy 
has also been studied [10, 12, 13]. In addition, it has been found that the accelerating 
cosmic expansion might be transient by performing a cosmographic evaluation [14] or 
using a bin method for q(z) [15]. 

More recently, using the SNIa (Constitution and Union2 [6, 7]), BAO distance ratio 
from the joint analysis of the Two-Degree Field Galaxy Redshift Survey (2dFGRS) and 
Sloan Digital Sky Survey (SDSS) data [8] and CMB shift parameter from Wilkinson 
Microwave Anisotropy (WMAP) [16, 17], Cardenas and Rivera [18] found that the tension 
between the low and high redshift data can be ameliorated effectively if the curvature Qk 
is incorporated as a free parameter, and both SNIa+BAO and SNIa+BAO+CMB favor 
a rapid variation in sign of q(z) between z ~ 0.5 and 0. 

In this paper, we plan to reanalyze the evolutionary behavior of q(z) with the curvature 
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Qk included as a free parameter in our discussion. The latest SNIa, BAO and CMB data 
will be used. Different from [18] where only BAO distance ratio is considered, we use 
four different BAO data sets to probe the systematics. These include the data at the low- 
redshift z = 0.106 [19] from 6-degree Field Galaxy Survey (6dFGS), at redshifts z = 0.35 
and 0.20 from the SDSS and the 2dFGRS [8], at redshifts z = 0.44, 0.6 and 0.73 from the 
baryon acoustic peak released by the WiggleZ Dark Energy Survey [20], and at z = 0.57 
from the Baryonic Oscillation Spectroscopic Survey (BOSS) [21]. Furthermore, for CMB 
data, we use the measurements of three derived quantities from WMAP7 [17]: the shift 
parameter R(z*) and the acoustic index l a (z*) at the recombination redshift z*, rather 
than only the shift parameter, since these three quantities can give unbiased information 
on dark energy parameters relative to the full CMB analysis as demonstrated in [22]. 

II. OBSERVATIONAL DATA 

For SNIa data, we use the Union2 data set released by the Supernova Cosmology 
Project (SCP) Collaboration [7], which consists of 557 data points. To constrain a theo- 
retical model from SNIa, we minimize the \ 2 value of distance modulus 



where fj,(z) = 5 log 10 [<i L (z) /Mpc] + 25 is the theoretical value of distance modulus, and 
fiobs is the corresponding observed one. di is the luminosity distance and it is defined as 



where the dimensionless Hubble parameter E{z) = H(z)/H , and Sk(x) is defined as x, 
sin(x) or sinh(x) for k — 0, +1, or — 1, respectively. For a dark energy model with the 
equation of state w(z), one has 

e 2 (z) = n m (i + z f + n k (i + zf + n de exp (3 J* l±^p.d/j , (3) 

where Q m and Qde are the dimensionless density parameters of matter and dark energy, 
and Vt de = 1 — Q m — Q k . Thus, for the CPL parametrization, w = w + Wiz/(1 + z), with 




(1) 
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w and Wi being two model parameters, we get 

E 2 (z) =n m (l + zf + fl fc (l + z) 2 + 0^(1 + z )3d-H»o+«;i) exp ( _ . (4) 

Since Ho is a nuisance parameter, we marginalize over it with a flat prior and then obtain 



2 r ^_^«? (E^i/g? ~ lnlO/5) 2 /In 

XSAT/alPJ - 72 1/ 2 I 



(5) 
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where a, = fi bs(zi) — 25 — 5 log 10 [iTo^L(^i)], and p denotes the fitting parameters in the 
model. 

For BAO data, we use the measurements from the 6dFGS (hereafter BA01) [19], 
the combination of SDSS and 2dFGRS [8] (hereafter BA02), the WiggleZ dark energy 
survey (hereafter BA03) [20] and the BOSS (hereafter BA04) [21]. For BAOl, the 
distance ratio d z = r s (zd)/D v (z) is measured at redshift z = 0.106 and the result is 
dp 106 = 0.336 ± 0.015 [19], where the effective distance is 



Dv(z) 

the drag redshift z^ is fitted as [23] 



[(1 + z) 2 H(z) 



1/3 



(6) 



with h = H /100, bi = 0.313(fi m /i 2 )-°- 419 [l + 0.607(fi m /i 2 ) - 674 ] and b 2 = 0.238(Q m h 2 ) - 223 , 
and the comoving sound horizon is 



where the sound speed is c s (z) = 1/yJ 3[1 + R b /(1 + z)\ with R b = 3fi 6 /(4 x 2.469 x 10~ 5 ), 
and Qb being the dimensionless baryon matter energy density. The constraint from BAOl 
can be obtained by minimizing 



_ (dp.106 ~ 0-336) 2 . . 

Xm ~ 0015^ 



From SDSS and 2dFGRS, Percival et al. [8] measured two distance ratios at redshifts 
z = 0.2 and z = 0.35 and obtained = 0.1905 ± 0.0061, d%% = 0.1097 ± 0.0036. Using 

2 

Xla= MiC^fadj) Adj, (10) 



one can obtain the result from BA02. Here Adi = di — d° bs and Cd z is the covariance 
matrix for two parameters c?o.2 and cio.35 given in [8]. 

BA03, containing three data points at z = 0.44, 0.6 and 0.73, is given by the WiggleZ 
dark energy survey [20]. The acoustic parameter rather than the distance ratio 

is measured in [20]. The constraint on model parameters can be obtained by using 

3 

x1 3 =J2 AA ^a( A ^ A ^ AA ^ (12) 

where the covariance matrix C^(Ai, Aj) is given in Table 2 in [20]. 

Recently, BOSS [21] released a BAO data A(0.57) = 0.444 ±0.014, which can be used 
to test theoretical model by minimizing 

, (AjpST) - 0.444) 2 

Xm ~ o3i4^ ' [ 6) 

Using Xbao = Xbi + Xb2 + Xb3 + Xba^ we can obtain the result from the combination 
of all BAO data sets. 

Except for the low redshift SNIa and BAO data, we also use the high redshift CMB 
information by implementing the Wilkinson microwave anisotropy probe 7 year (WMAP7) 
data. Since, when the full WMAP7 data are applied, some more parameters depending 
on inflationary models need to be added, the ability to constrain dark energy models will 
be limited. Thus, we use three derived quantities [17]: the shift parameter R(z*) and the 
acoustic index Ia{z*) at the recombination redshift z*. 

fl(2 "^ftK^fi|))^ 725±0 - 018 - (14) 

''"•> = (ilttl-r m09±076 - (15 » 

where the redshift z* is given in [24] 

z* = 1048[1 + 0.00124(fi fe )-°' 738 ][l + gi (n m h 2 ) 92 ] = 1090.04 ± 0.93 (16) 
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with 



_ 0.0783(^ 2 )- - 238 
91 l + 39.5(fi 6 /i 2 ) - 763 ' [ ] 



0.560 

92 



l + 21.1(fi 6 /i 2 ) 1 - 81 v ' 

Thus, the Xcmb can be expressed as 

Xc*M_b(p) = ^ij^l^^CcMB^ii Xj)AXj, (19) 

where the three parameters = [R(z*), l^(z*), z*], Axj = Xi — x° bs and covariance matrix 
CcMBi^i, Xj) is taken from Table 10 in [17]. It has been found that these CMB quantities 
can give similar constraints on dark energy parameters compared with the full CMB power 
spectrum [22]. 



III. RESULTS 



We use the Monte-Carlo Markov Chain (MCMC) method to explore the model param- 
eter space p. The MCMC method randomly chooses values of parameters p, calculate % 2 
and determines whether to accept or reject the set of parameters p using the Metropolis- 
Hastings algorithm method [25]. 

To probe the properties of dark energy and the cosmic expansion history, we analyze the 
constraints on the CPL parametrization and then reconstruct the evolutionary behavior of 
the decelerating parameter q(z) = —da/ a 2 . The q(z) gives information on the expansion 
speed and acceleration since its sign tells us whether the expansion is accelerating or 
decelerating, and its first derivative shows whether the acceleration is slowing down or 
speeding up. 



A. spatially flat case 

The constraints on model parameters Q m , Wq and W\ are given in Tab. (I) and Fig. (1), 
and the evolutionary behaviors of q(z) are shown in Fig. (2). It is easy to see that the 
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Union2 SNIa alone only has a weak constraint on model parameters. The best fit line of 
q(z) shows that a slowing down of expansion acceleration, whose phase transition occurs at 
z ~ 0.3, is favored by SNIa. Adding BAO data into the analysis, one can find that BAOl 
and BA02 have an almost negligible effect on results. However, when BA04 is included, 
the model parameters, especially w±, are tightened apparently, but a slowing down of the 
expansion acceleration is still preferred. Once BA03 is considered, a tighter constraint is 
obtained and, different from the cases of BAOl, BA02 and BA04, BA03+SNIa favors 
that the cosmic acceleration is speeding up. Thus, different BAO data sets seem to give 
different results, which implies that there exists a tension between them. Combining all 
BAO data and SNIa data (SNIa+BA01+BA02+BA03+BA04), we find that the phase 
transition from speeding up to slowing down remains to be favored, but the transition 
redshift is very close to zero. 

When the WMAP7 CMB data is included, the slowing down phenomenon of the expan- 
sion acceleration, obtained in BAOl+SNIa, BA02+SNIa, BA04+SNIa and BAO+SNIa, 
disappears in agreement with what was obtained in Ref [11, 12, 18] where the CMB shift 
parameter and BAO distance ratio between redshift 0.2 and 0.35 are used. However, 
BA03+SNIa and BA03+CMB+SNIa give a very consistent result. So, the tension be- 
tween the low redshift BAO+SNIa data and the high redshift CMB data may not always 
exist. In addition, BA04+CMB+SNIa seem to show a slight difference on constraining 
model parameters and it allows larger Q m and Wi, and a more negative Wq. This difference 
can also be seen from the evolutionary curve of q(z) in the middle panel of Fig. (2). 

B. spatially curved case 

In Ref. [18], it has been found that the spatial curvature has a significant effect on 
reconstructing the evolution of q(z) and alleviating the tension between low redshift data 
and high redshift one since both SNIa+CMB shift parameter+BAO distance ratio, and 
SNIa+BAO distance ratio favor a slowing down acceleration when the spatial curvature 
is considered, which is just opposite to the flat case. Thus, here, we also study the effect 
of spatial curvature. The results are shown in Tab. (II) and Fig. (3). The corresponding 
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q(z) evolutionary curves are given in Fig. (4). We find that SNIa and SNIa+BAO (includ- 
ing SNIa+BAOl, SNIa+BA02, SNIa+BA03, SNIa+BA04 and SNIa+BAO) give very 
consistent constraints on model parameters. This differs from the flat case where BA03 
leads to a different result. An open universe is favored slightly although Qk = is still 
allowed at the la confidence level. Therefor, the tension between different BAO data sets 
obtained in flat case disappears when the spatial curvature is taken into account. The left 
and right panels of Fig. (4) show that the cosmic acceleration has entered a slowing down 
era and q > seems to be favored, which means that the accelerating cosmic expansion 
may be a transient phenomenon and the universe may have re-entered the decelerating 
phase. In addition, we find that a very large Wq, wo ~ —0.5, is given by SNIa and BAO 
data. 

When the WMAP7 CMB data is further added into our analysis, the best fit results 
change drastically. In this case, is very close to zero and thus a flat universe is 
supported strongly. Except for the SNIa+BA04+CMB case, which favors a phantom- 
like dark energy at the present with a crossing of —1 line in the near past, observations 
prefer a Lambda cold dark matter model since wq is very close to —1 and w\ is very 
small. The middle and right panels of Fig. (4) show that a speeding up expansion is 
supported when CMB is included, which is the same as the flat case, but is different 
from what was obtained in Ref. [18] where only the CMB shift parameter is considered. 
Thus, comparing the results in spatially flat and curved cases, we find that the inclusion 
spatial curvature leads to somewhat more serious tension between SNIa+BAO and CMB, 
although it alleviates it between different BAO data sets. 

IV. CONCLUSION 

In this paper, using the MCMC method, we reconstruct the evolutionary behavior of 
the decelerating parameter q(z) from the latest observational data including the Union2 
SNIa, BAO, and CMB data from WMAP7. For BAO data, four different kinds of data ob- 
tained from the 6dFGS, the combination of SDSS and 2dFGRS, the WiggleZ dark energy 
survey and the BOSS are used. In our analysis, the CPL parametrization for the EOS 
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TABLE I: The marginalized la constraints on flat CPL model by different observational data. 



Data set 






w 




SN 


541.4544 


416n +al054 




. 9981+6,9081 
°- zzoi -19.0726 


SN+BAOl 


541.4586 


4162+ ' 1105 

U.41OZ_ 2901 




c: 4 r;oo+6.9106 
0.405Z_2o.5981 


SN+BA02 


542.1974 


4246+ ' 0925 
UA/w -0 .2780 


-0.8600i ;^ 


c oi oi +7.1387 
o.oioi„ 18 5244 


SN+BA03 


542.6020 


2Q23+ a0550 


-1-10191°:^ 


-0.2590l 2 ;^? 


SN+BA04 


541.9226 


0-3580+S 


-1.01021S 


-1.9484±|; 4 ^1 


SN+BAO 


549.3644 


0.31681°;^! 


-0.9689l°|f 7 ° 


-l iQin +2 - 0452 

i.iyiu_ 2 7636 


SN+BAOl+CMB 


543.2798 


0.2789^ 


-i.oi78±g : gg 


o.ooiii?S 


SN+BA02+CMB 


544.2428 


0.2772±8;g^ 


-i.o39o±8:2gg 


0.11451?;™ 7 


SN+BA03+CMB 


543.1966 




-1 0323+ ' 2477 

1.UO^O_q 2094 


u - uu ' z -1.4231 


SN+BA04+CMB 


548.4828 


n ^901 +0-0294 

u.ozui_ 0487 


-1 41 42+0-6588 

l.*l^Z_Q 14g7 


1.49041™ 


SN+BAO+CMB 


554.8478 


2939+ - 0269 
u.zyoy_ 0231 


-1 0464+ - 2848 
LU ™-0 .3129 


-0.21061^ 



of dark energy is considered. For a spatially flat universe, we find that there is tension 
between different BAO data sets since BAOf +SNIa, BA02+SNIa and BA04+SNIa sup- 
port a slowing down of the acceleration of the cosmic expansion, while BA03+SNIa does 
not. When the WMAP7 CMB data (R(z*), Ia{z*) and z*) is added into our discussion, 
the slowing down phenomenon disappears and rather a speeding up is favored. Thus, 
there is also tension between BAO+SNIa and CMB, except for the case of BA03 since 
BA03+SNIa and CMB+BA03+SNIa give a very consistent constraint on model param- 
eters. By incorporating the spatial curvature as a free parameter, we find that SNIa and 
SNIa+BAO seem to support an open universe, the accelerating cosmic expansion may 
be a transient phenomenon, and the tension between different BAO datasets is allevi- 
ated effectively. However, when the CMB data from WMAP7 is included, observations 
favor strongly a Lambda cold dark matter model and a spatially flat universe. Mean- 
while, observations with CMB included prefer that the cosmic acceleration is speeding 
up, which is consistent with the flat case and what was obtained in Refs. [10, 11], but 
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FIG. 1: The la and 2a contours in wq — Cl m (a) and wq — w\ (b) planes for the flat CPL 
model. The dashed, solid and thicken solid lines shows the results from SNIa, SNIa+BAO and 
SNIa+BAO+CMB, respectively. 



is different from what was obtained in [18] where it was found that in a spatially curved 
universe, the SNIa+BAO distance ratio+CMB shift parameter gives a consistent result 
with SNIa+BAO distance ratio and both of them favor a slowing down of the cosmic ac- 
celeration. Comparing our results obtained in the flat and curved cases, one can see that, 
although the spatial curvature can alleviate effectively the tension between different BAO 
data, it leads to somewhat more serious tension between SNIa+BAO and CMB. After all, 
in a flat universe, SNIa+BA03 and SNIa+BA03+CMB give a very consistent constraint. 
This result is in sharp contrast to that given in [18] where only a CMB shift parameter is 
considered. Since R(z*), Ia{z*) and z* give similar constraints on dark energy parameters 
compared with the full CMB power spectrum as is shown in [22], we think that our results 
may be more reliable than that reached in [18]. 
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(a) (b) (c) 




FIG. 2: The best fit evolutionary curves of the deceleration parameter for the flat CPL model. 
In (a), the solid, dashed, dot-dashed and dotted lines represent the results obtained from 
SNIa+BAOl, SNIa+BA02, SNIa+BA03 and SNIa+BA04, respectively. In (b), the solid, 
dashed, dot-dashed and dotted lines represent the results obtained from SNIa+BAOl+CMB, 
SNIa+BA02+CMB, SNIa+BA03+CMB and SNIa+BA04+CMB, respectively. In (c), the 
solid, dashed and dotted lines show the results from SNIa, SNIa+BAO and SNIa+BAO+CMB, 
respectively. 
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